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Abstract 

We formulate stochastic partial differential equations on Riemannian manifolds, mov- 
ing surfaces, general evolving Riemannian manifolds (with appropriate assumptions) and 
Riemannian manifolds with random metrics, in the variational setting of the analysis 
to stochastic partial differential equations. Considering mainly linear stochastic partial 
differential equations, we establish various existence and uniqueness theorems. 
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1 Introduction 

Stochastic partial differential equations (SPDE) are becoming increasingly popular in the math- 
ematical modelling literature. Analogous to the difference between ordinary differential equa- 
tions (ODEs) and partial differential equations (PDEs), it seems that in some cases, stochastic 
differential equations are not as accurate as describing physical phenomena as SPDEs are. It 
is because of this, and the want of generalising ltd diffusions to infinite dimensions for applica- 
tions to problems in physics, biology and optimal control that the theory of SPDEs has grown 
exponentially in the past four decades. 

However, there seems to be a distinct lack of mathematical theory for SPDEs on moving 
surfaces, at odds with the deterministic counterpart. Indeed, a survey into the mathematical 
literature for SPDEs on moving surfaces produces no results. The use of such objects is wide- 



spread in the applied literature (Meinhardt 1982, 1999 ; Neilson et al. 2010 amongst others) 
and indeed the paper by Neilson et al. 2010 along with the suggestion of Professor Charles 



Elliott prompted this study into the objects. 

If we go one step back and ask for SPDEs on (Riemannian) manifolds, instead of moving 
surfaces, we find three papers |Gyongy| [l993| (1997) and |Funaki| (1992) . The last paper considers 
SPDEs whose solution is a function / : S — >■ M. where S is the unit disc and M. is the 
manifold. Although such objects are prevalent in mathematical physics (Funaki 1992 and 



references within) we are only interested in SPDEs whose solution is a real valued function 
g : M. — > K. Indeed, the only theory for SPDEs on manifolds with real-valued functions as 



solutions is given in Gyongy 1993, 1997 



There are three main approaches to analysing SPDEs, namely the "martingale approach" (cf 
Walsh] [1986] ), the ' 'semigroup (or mild solution)" approach (cf |Da Prato and Zabczykj |1992| ) 
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michael.r.scott@warwick.ac.uk. The research for M.S was funded by EPSRC as part of the MASDOC DTC 
with grant reference number EP/H023364/1. 
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and the "variational approach" (cf Rozovskii 1990 , Prevot and Rockner 2007|). The ap proach 



of SPDEs on a differentiable manifold in Gyongy 1993 is that of Da Prato and Zabczyk 1992 
namely the semigroup approach. 

There is no mathematical literature for the variational approach to SPDEs on Riemannian 
manifolds and for this reason, we adopt this approach in this paper. Here we pose and give ex- 
istence and uniqueness results for SPDEs on Riemannian manifolds, SPDEs on moving surfaces 
and finally SPDEs on evolving Riemannian manifolds, which allows us to look at Riemannian 
manifolds with random metrics. This paper is organised in the following way: 

In chapter [2] we proceed to define stochastic partial differential equations in the general 



setting. This will be an abstract setting and where we mainly follow the monograph of Prevot 



and Rockner |2007| . After giving notation and elementary definitions, we briefly look at the 



abstract definition of what a SPDE is in terms of the variational approach, giving an existence 
and uniqueness result, concluding by giving an example. 

In chapter [3] we formulate what it means to have an SPDE on a Riemannian manifold, 
M.. We give a self-contained (presenting results without proof) introduction to Riemannian 
geometry which sets up all the necessary theory to define differential operators for smooth 
functions / : M. — > R. Following this, we define the Sobolev spaces needed and prove the 
Poincare inequality which is needed for a later example. Having set all the preliminary theory, 
we define what it means to have a SPDE on a Riemannian manifold and consider two specific 
examples; proving an existence and uniqueness result in each case. For the examples we con- 
sider the stochastic heat equation whilst the second example is the non-degenerate stochastic 
heat equation, where the Laplace-Beltrami operator is replaced with the p— Laplace-Beltrami 
operator, for p > 2. 

In chapter [4] we study SPDEs on moving hypersurfaces. Firstly, we define what we mean 
by a hypersurface giving all the necessary theory. Following this we formulate a deterministic 
PDE on a moving surface as a consequence of conservation law which allows us to consider 
the stochastic analogue (which includes choosing the noise) of this object. This turns out to 
be the stochastic heat equation on a general evolving hypersurface Ai(t). We always assume 
that A4(t) is compact, connected, without boundary and oriented for all t G [0,T], with points 
evolving with normal velocity only. Penultimately we consider the concrete example of when 
M.(t) is the 5* n_1 sphere evolving according to "mean curvature flow" and we finally consider 
the nonlinear stochastic heat equation on a general moving surface, where points on the surface 
evolve with normal velocity only, noting that the nonlinearity is not in any of the derivatives. 

In chapter [5] we change how we think about a manifold evolving. Instead of thinking of a 
one-parameter family of manifolds M(t), t G [0, T] we think of one manifold Ai with a one- 
parameter family of metrics g(-,t), t G [0, T\. As given in the discussion section of chapter |5j we 
will see that under specific technical assumptions, the equations that live on M. are equivalent 
to the equations that live on Ai(t). We also see that this change of view enables a more natural 
noise to be chosen, as supposed to the one chosen in chapter [4} Following the discussion, we give 
an existence and uniqueness theory for a general parabolic SPDE, with minimal assumptions 
for which the approach works. Following this, we consider a random perturbation of a given 
initial metric, which we will refer to as a "random metric". 

Chapter[6]is the final chapter, detailing possible extensions to this paper for further research. 
We detail the mathematical challenges needed to be overcome in order to solve the problems 
outlined. 

Thanks go to C.M.E and M.H for supervising M.R.S during this project. 
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Stochastic partial differential equations: The general 
setting 



2.1 Notation and definitions 



We adopt the notation and give definitions as in Prevot and Rockner 2007 . Throughout 
this paper we fix T G (0, oo) and a probability space (Q, J 7 , P) with filtration T% that satisfies 
the usual conditions, i.e it is right continuous and J-q contains all the P— null sets. For X a 
(separable) Banach space, we denote by B(X) the Borel a- algebra. Unless otherwise stated, all 
measures will be Borel measures. 

Let if be a separable Hilbert space with inner product (•, •)# and induced norm || • \\ H . 
Suppose V is a Banach space with V C H continuously and densely. By this we mean there 
exists C > such that \\u\\h < f° r every u G V and that given u G H there exists a 

sequence Uk G V such that \\v,k — u\\h — > as k — > oo. For the dual of V, denoted V* := {I : 
V — > R 1 1 linear and bounded} we have that H* C V* continuously and densely and identifying 
H and H* via the Riesz isomorphism we have 

V C H C V*. 

Such a triple is called a Gelfand triple. We denote the pairing between V* and V as (•, •) and 
note that for h G H and v G V we have (h, v) = (h, v)h- 

We denote by L(X, Y) all the linear maps from X to Y . When X = Y we write L(X) 
instead of L(X,X). 

If X and Y are separable Hilbert spaces and {ej}"^ is an orthonormal basis of X then 
T G Y) is called Hilbert-Schmidt if 

ll T lli 2 (x,y) : = J] ( Tei ' Te ^ < 00 ( 2>1 ) 



/el 



and is called finite-trace if 



tr (T) := X (Te^Ci) < oo. 



We denote the linear space of all Hilbert-Schmidt operators from X to Y by L,2(X,Y) and 



equip this space with the norm defined in (2.1) 



Fix U a separable Hilbert space, T G (0, oo) and Q G £(£/) such that Q is non-negative 
definite, symmetric with finite trace (which implies that Q has non-negative eigenvalues). 

Definition 2.1. A U— valued stochastic process W{t), t G [0, T], on a probability space (ft, J 7 , P) 
is called a (standard) Q— Wiener process if 

1. W(0) = 0; 

2. W has P— a.s continuous trajectories; 

3. The increments of W are independent. That is, the random variables 

W(h), W(t 2 ) - W(h), ■■■ , W(t n ) - W{t n ^) 
are independent for all < t\ < ■ ■ ■ < t n < T, n G N; 
4- The increments have the following Gaussian laws 

Fo(W(t) -W(s))- 1 =N(0, (t-s)Q) for every < s < t < T. 
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Note that the definition of the stochastic integral can be generalised to the case of cylindrical 
Wiener processes, where the covariance operator need not have finite trace. The reader is 
directed to Prevot and Rockner 2007 for a more general discussion. 



2.2 Abstract theory of stochastic partial differential equations 

In the following we will fix U a separable Hilbert space and let Q = I. Let W be the resulting 
cylindrical Wiener process. It is this object that will be the mathematical model for "noise" in 
the SPDEs. 



We will follow Prevot and Rockner 2007 chapter 4 for the formulation, statements of the 
existence and uniqueness theorem and their consequent proofs. 

Let H be a. fixed separable Hilbert space with inner product (•, and denote by H* its 



dual. Let V be a Banach space such that V C H continuously and densely as in section 2.1 



Consider the Gelfand triple V C H C V* as discussed in section 2.1 Here, B(V) is generated 
by V* and 13(H) by H*. 

We wish to study stochastic differential equations on H of the type 

dX(t) = A(t, X(t)) dt + B(t, X(t)) dW(t) 

x(o) = x {2 - 2) 

where X is a given stochastic process. 



We will refer to such equations (2.2) as stochastic partial differential equations (SPDE) for 
when A is a differential operator. 

The important point to realise is that an SPDE is an infinite dimensional object. It is quite 
useful to think of such objects as "PDE + noise". Indeed, even though A and B take values in 
V* and L 2 (U,H) respectively, the solution X will, however, take values in H again. For when 
V and H are function spaces and A is a differential operator this means that the solution is 
function valued, which is perhaps a difficult concept to comprehend at first. 

We proceed to give the precise conditions on A and B that will be considered through the 
paper. 

Fix T G (0, oo) and let (Q, J 7 , P) be a complete probability space with normal filtration J- t , 
t G [0,T]. We assume that W is a cylindrical Q- Wiener process with respect to J^, t G [0,T], 
taking values in U and with Q = I. 
Let 

A: [0,T] x V Xfi — ► V* 

B : [0,T] x V x Q — ► L 2 (U,H) 

be progressively measurable. By this we mean that for every t G [0,T] the maps A and B 
restricted to [0, t] x V x Q are B[0, t] £g> B(V) <S> J-i-rneasurable. When we write A(t, v) we mean 
the map u — > A(t,v,u) and analogously for B(t,v). 

Assumption 2.2. The following hypotheses will be on A and B throughout the paper. 
(HI) (Hemicontinuity) For all u, v, x G V, u G f2 and t G [0,T] the map 

K 3 A H> (A(t, u + \v, u),x) 

is continuous. 
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(H2) (Weak Monotonicity) There exists cGf such that for every u, v G V 

2{A(-,u)-A(-,v),u-v) + \\B(-,u)-B 
< c\\u-v\\ 2 H on [0,T] x fl 



v )\\l 2 (u,h) 



(H3) (Coercivity) There exists a G (l,oo), c\ G R, c 2 G (0, oo) and an (J 7 *)- adapted process 
f G ^([0, T] x tt, dt g> P) snca too* /or every n G V, t G [0, T] 



2(A(t,n),n) 



W^)lli aW H)<ci|Hllr-c 2 |Hlv + /(*) 



on f2. 



(H4) (Boundedness) There exists c 3 G [0, oo) and an [Tt)-adapted process g G L«- 1 ([0,T] x 
Cl, dt g) P) snca taat /or every n G V, t G [0, T] 



\\A(t,v)\\ v * <g(t)+c 3 \\v 



|a-l 



where a is the same as in H3. 

These hypotheses appear to be quite abstract and on the face of it, and so we give some 
intuition as to why they are needed. 

One can see that H3 and H4 really come from the deterministic case of the variational 
approach to PDE (Evans 1998| ). Note also that in the case of A being non-linear, H2 is also 
common in deterministic PDE theory. Indeed, the method of Minty and Browder (Renardy 



2004] ) uses the monotonicity of A to identify the weak limit of A(uk) as A(u) (here Uk is some 



Galerkin approximation to the solution u). Furthermore, as in the case of Minty and Browder, 
continuity of u i— > A(-,u) is used and so HI is a natural generalisation of this. 

The reader should observe that as soon as A is linear on V, HI is immediately satisfied by 
the definition of the pairing between V and V*. To see this let t G [0, T] and u,v,x G V and 
u> G Q. Then 

(A(t, u + An, u),x) = (A(t, u, uj),x) + X(A(t, v, oj),x) 

and so R 3 A i— > (A(t, u + An, u), x) is clearly continuous. 

Later we will give examples of A and B and of the spaces V, H and V* but first we proceed 



to define exactly what we mean by "solution" to (2.2), as taken from Prevot and Rockner 2007 
page 73. 

Definition 2.3. A continuous H-valued {Ft) -adapted process X(t), t G [0, T], is called a so- 



lution of ( |2T2| ) ; if for its dt <g> F- equivalence class X we have X G L a ([0,T] x Q, dt ® P; V) n 
L 2 ([0,T},n,dt<S)F;H) with a as in H3 andF-a.s 

X(t)=X{0)+ [ A{s,X{s))ds+ [ B(s,X(s))dW(s), t G [0,T], 
Jo Jo 

where X is any V -valued progressively measurable dt ® F-version of X- 

For the technical details of the construction of X, the reader is directed to exercise 4.2.3 of 



Prevot and Rockner 2007 page 74 



The following is the main existence result, which was originally proven in |Krylov and R(> 
Instead of giving the proof in its entirety, we outline the ideas and refer the 



zovskii 1979 



reader to the relevant pages of Prevot and Rockner 2007 



Theorem 2.4. Let A and B satisfy assumption 2.2 and suppose X G L (Q, J^o^F; H). Then 



there exists a unique solution X to (|2.2|) in the sense of definition 2.3\ Moreover, 

\X( 



E 



sup 

te[o,T] 



< oo. 
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2.3 An example 

In the following we give a concrete example of an SPDE on A C lR n where A is open and the 
boundary <9A is sufficiently smooth for the required Soblev embeddings. The following example 



is taken from Prevot and Rockner 2007 , and the reader is referred to this text (pp. 59-74) for 
further examples. 

Let A = A := Y17=i ih 1 ^ e ^ ne Laplacian, with domain 



Cg° (A) 



{u G C°°(A) : supp(w) is compact}, 



where supp(w) = {x G A : u(x) ^ 0}. Recall the Sobolev space (Adams 2003 ) -f^ 1 (A) := {u : 
A — > R : u G L 2 (A), |Vm| G L 2 (A)} where Vw exists in the weak sense, and that Hq(A) is 
defined as the closure of C^°(A) in the norm 



ml hi 



\u 



|2 , 

II 2 + 



IVttl 



|2 

II 2 - 



|2 

II 2 - 



To save on typesetting, we will abuse notation and write ||Vu||£ 2 for || |Vm| 

It is well known (Evans 1998| ) that A has a unique extension from Qj°(A) onto Hq(A). 
Thus, define V = Hq(A) and observe that V C L 2 (A) continuously and densely (Evans |1998|, 
Sobolev embedding). Define H := L 2 (A) and identifying H with its dual H* we will consider 
the Gelfand triple V C H C V*, or more concretely i^o(A) C L 2 (A) C H~ x , recalling the 

:=(itf(A))*. 



Evans 1998 that H 



notation m 

So we have defined the operator A and the associated Gelfand triple. For the noise, we fix 
U some abstract separable Hilbert space, and ask for some Hilbert-Schmidt map from U to H. 
It is not important what U is, for if i : U — > H is Hilbert-Schmidt and time independent, the 
noise idW interpreted as the stochastic integral J^idW(s) which lies in H. We have 

Proposition 2.5. Let U and H be fixed separable Hilbert spaces. Then there exists i : U — >■ H 
which is Hilbert-Schmidt. 

Proof. Let {ej} ig N, {/j}j 6 N be orthonormal bases of U and H respectively, which exist as U and 
H are both separable. Define 



Then % : U — > H is Hilbert-Schmidt since 



(*(ei),i(ej))n 



which is summable. 



□ 



From this, we will fix U some abstract separable Hilbert space^j] and % : U — > L 2 (A) Hilbert- 
Schmidt as constructed in proposition 2.5 and consider the SPDE on A which we will call the 
"Stochastic Heat Equation" 



dX{t) = AX{t) dt + idW{t) 
X(0) = X 



(2-3) 



1 Indeed from proposition 2.5 one can take U = H = L 2 (A) 



6 



where X G L 2 (fl, Fq, P; H) is given. Note here (0, J 7 , P) is a complete probability space and 
the cylindrical Wiener process W(t), t G [0,T], is with respect to a normal filtration (J 7 *). 

Note here that A and i do not depend on the probability space (Q, J 7 , P) and so trivially i 
is predictable and since it is Hilbert-Schmidt, the stochastic integral j^idW{s) is well defined. 

We have the following 



Proposition 2.6. Suppose that X G L 2 (Q, Fq, P; H). Then the stochastic heat equation (2.3) 
has a unique solution in the sense of definition 2.3[ Furthermore, 

E[ sup < oo. 

*e[o,T] 



Proof. From theorem 2.4 it suffices to show that A := A and B := i satisfy HI to H4 of 



assumption 2.2 



1. Since A is linear we see that HI is satisfied. 

2. To see H2, observe that as i is independent of the solution X we have that for any u,v G V 
\\B(-,u) — B(-,v)\\l 2 (u : h) = 0. Also, by definition of V, there exists Uk,Vk G Cq°(A) such 
that Uk — > u and v & — >■ v in V. Hence 

2(Au — Ai>, u — v) = 2 lim (Am^ — Au*., Mfc — v/c)i? 
= -2 lim ||V(w fc - 
= -2||V(«-t;)||3! r 

2 „ 1.2 

where C p is the Poincare constant from the Poincare inequality (Adams 2003] ) which 
says that there exists C v > such that for every u G Hq{K) 

\\u\\ h < CpllVwIln- 
Thus H2 is satisfied with c = —2/C 2 . 

3. To see H3, using the same argument as above for v G V 

2(Av,v) = -2\\Vv\\ 2 H = 2\\v\\ 2 H - 2\\v\\ 2 v 

since \\v\\y = |H|g + ||Vf \\%. Since i is Hilbert-Schmidt there exists k > such that 
\\i\\L 2 (u,H) < k, hence 

2(Av, v) + IHH!^) < 2||^||1, - 2\\v\\ 2 v + k 2 . 

Noting that k 2 is (J-t)-adapted and is in -^ 1 ([0, T] x Q, dt <8> P), we see that H3 is satisfied 
with a = C\ = c 2 = 2 and f(t) = k 2 . 

4. Finally, for H4, if u,v G C^(A) then 

\{Au,v}\ = \(Au,v) H \ = \(Vu,Vv)h\ < ||Vw||if||Vu||H < 

which implies that ||Au||y* < for every v G V by a density argument, and so H4 is 
satisfied with c 3 = 1 and g(t) =0. 



Now applying theorem 2.4 we see that (2.3) has a unique solution. □ 

Remark 2.7. In item 2 above, we have that 2(Au — Av,u — v ) = — 2||V(« — v)\\^j < and so 
we could have taken c = for H2. Thus, there is no need to use the Poincare inequality. This 
point will be important later. 
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3 Stochastic partial differential equations on Rieman- 
nian manifolds 



3.1 A brief introduction to Riemannian manifolds 



In order to define what we mean by SPDEs on Riemannian manifolds, we must first have a 
working knowledge of the theory of Riemannian manifolds. This is referred to as Riemannian 
geometry in the literature. 

There are many introductory texts to Riemannian manifolds such as Lee 1997, 2003 and 



Hebey 1996 chapter 1. For a more advanced text in general differential geometry the reader 



is directed to Spivak 1999 . We first introduce smooth manifolds as in Lee 2003 



Definition 3.1. We say Ai is a smooth manifold of dimension n if Ai is a set and we are 

given a collection {U a } a of subsets of Ai together with an injective map : U a — > R n for each 
a such that the following hold: 

1. For each a, the set <f(U a ) is an open subset ofR n ; 

2. For each a, (3 the sets (p a (U a H Up) and ipp(U a fl Up) are open in W 1 ; 

3. Whenever U a nU/3 ^ the map ^o^ 1 : (pp(U a nUp) — > <£ a {U a nUp) is a diffeomorphism; 
4- Countably many of the sets U a cover Ai; 

5. For p 7^ q where p, q G Ai either there exists U a with p,q G U a or there exists disjoint 
U a , Up such that p & U a and q G Up. 

We say that each (U a ,ip a ) is a smooth chart; that is U a C Ai is open and <p a : U a — > <p a (U a ) 
is a homeomorphism. 

We will need some notion of smoothness for functions / : Ai — > R. The notion of smooth- 
ness for such / is inherited from the notion of smoothness of functions g : R n — > R. Precisely: 

Definition 3.2. Let Ai be a smooth manifold. We say f : Ai — > R is smooth if for every 
p G Ai there exists a smooth chart (U, tp) for Ai whose domain contains p and such that 
f o if -1 : ip{U) — > R is smooth on the open subset U := ip(U) C R n . 
The set of all such functions will be denoted by C°°(Ai). 

An important observation is that Ai is not a vector space in general. For example, if one 
takes Ai := S 11 ' 1 := {x G lR n : \\x\\ = 1} then if x,y G Ai then + = 2 and so x + y ^ Ai. 
However, to each point p G Ai there is an associated vector space structure. This is referred 
to as the tangent space. 

Definition 3.3. Let Ai be a smooth manifold and let p G Ai. A linear map X : C°°(Ai) — > R 
is called a derivation at p if X(fg) = f(p)X(g) + g(p)X(f) for every f,g£ C°°(Ai). The set 
of all such derivations at p is called the tangent space at p and will be denoted by T p Ai . 



Observe that T p Ai is indeed a vector space. Further, it is shown in Lee 2003 page 69 that 



T p Ai is an n-dimensional vector with basis 





where the x l are local coordinates. 

Related to the tangent space is the so called tangent bundle. 
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Definition 3.4. We define the tangent bundle, denoted TM as 

TM := |J T P M, 

noting that this is a disjoint union. 

This now allows us to define the manifold analogue of a vector field. 

Definition 3.5. A vector field Y : M — > TM, usually written p Y p is such that Y p G T P M 
for each p G M. 

The set of all such vector fields will be denoted by C°°(M,TM) . 
Remark 3.6. Indeed, since T P M is a vector space, one has that 

d 



dx i 



where Y l : U — > "R (1 < i < n) are called the component functions of Y in the given chart 
(U,<p). 

With these constructions, it is natural to define a metric on T P M. 

Definition 3.7. Let g p : T p M x T P M — > M be symmetric and positive definite at each p G M, 
which means that g(u,v) = g(v,u) for every u,v G T P M and g(u,u) > for all u G T P M. 
Then g is called a metric on T p M. 

Remark 3.8. Since g is symmetric and positive definite, this leads to a positive definite and 
symmetric matrix (g^) G IR nxn defined via 



9ij ■= 9(di,dj) l<i,j<n 



where = . We refer to g^ as the components of the metric g. 



We now have all the theory to define a Riemannian manifold. 

Definition 3.9. A Riemannian manifold is a pair (M,g) where M is a smooth manifold and 
g is a metric. 

Remark 3.10. One can show using partitions of unity that given a smooth manifold M there 
always exists a metric g on M. The arguments are omitted. 

We will now write M for (M,g) and only consider Riemannian manifolds without boundary. 

In order to define SPDEs on M we will need to define differential operators on M. Further, 
to specify function spaces, we need some notion of integration on M. This will ultimately, in 
section |3T2| enable us to define Sobolev spaces on M. 



A step towards looking at differential operators on M is the notion of connection (Lee 1997 
page 49). 

Definition 3.11. A connection on M is a bilinear map 

C°°{M,TM) x C°°(M,TM) ^C°°{M,TM) 

(X,Y)^V X Y 

such that 
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1. VxY is linear over C°°(Ai) in X, that is 

V fXl+g x 2 Y = fV Xl Y + gV X2 Y for every f,ge C°°(M); 

2. VxY is linear over K in Y , that is 

Vx(aY! + bY 2 ) = aVxYx + bV x Y 2 for every a, b G K; 

3. V satisfies the following product rule 

Vjr(/y) = /V*y + (X/)y /or every / G C°°(A<) X, Y G C°°(.M,TA<). 



Analogous to remark 3.8 letting X = di and Y = dj we have 
Definition 3.12. 

V di dj = T™d m 

and we refer to T™ as the Christoffel symbol of the connection V. 

We will be considering a special type of connection on Ai; the Levi-Cevita connection. 

Theorem 3.13 (Fundamental theorem of Riemannian geometry). Let Ai be a Riemannian 
manifold. Then there exists a unique connection V on A4 that is compatible with the metric g 
and is torsion free. By this we mean that for every X,Y, Z G C°°(A4, TAi) 

Vxg{Y, Z) = (compatible with the metric) 
and VxY — VyX = [X, Y\ (torsion free). 

Such connection is called the Levi-Cevita connection. 



Proof. The reader is directed to Lee 1997 page 68 for the proof. □ 

We now define some differential operators that will be used. We define the gradient of a 
function u : M. — V R, denoted V«, as having representation in local coordinates 

(Vti)i = diU, 



noting that |Vu| 2 = g^diudjU (Hebey 1996 page 10) in local coordinates. We define the 



Laplace-Beltrami operator, Am, of a function u : M. — > K as 

n , 

A M u := -yr= { d m (V\i\g mk d k u) (3.1) 



k,m=l 



in local coordinates, where \g\ = det(py) and o ,J is the (i,j) th element of (oy) 1 , the inverse of 
(9ij)- 

Finally, for integration, one defines the Riemannian volume element 



df(flf) := a/M dx 

where dx is the Lebesgue volume element of lR n . 

The reader should note that we have not mentioned all the aspects of Riemannian geometry 
and in particular we have not mentioned curvature. We will not mention the various types of 



curvature one can define on Ai but refer the interested reader to Lee 1997 



In the next section we will introduce Sobolev spaces on Ai and give the precise assumptions 
that we will employ on Ai. This will setup the theory needed to define SPDEs on Ai. 
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3.2 Formulation of a stochastic partial differential equation on a 
Riemannian manifold 



The abstract theory of chapter 2 and the preceeding theory of Riemannian manifolds will now 
allow us to consider SPDEs on Ai. Analogously to section |2.3[ in order to define what we 



mean by a SPDE on a Riemannian manifold Ai, one needs to identify the differential operators 
acting on real-valued functions defined on Ai and the appropriate Gelfand triple. 

Essentially the only hard work one needs to worry about is whether or not the Sobolev 
embeddings that hold on an open subset of lR n (with a sufficiently smooth boundary), also hold 
on Ai. 

The topic of Sobolev embeddings on Ai is far from trivial. It turns out that many of the 
Sobolev embeddings that hold on IR n are simply false on a general Riemannian manifold. Two 
useful texts for Sobolev spaces on Riemannian manifolds are Hebey 1996, 2000 , but the work 



in this area arguably dates back to Aubin 1976 



For technical reasons, we employ 

Assumption 3.14. Ai is a compact Riemannian manifold of dimension 1 < n < oo 
which is connected, oriented and without boundary. 



Such an example of Ai is S n 1 
following. 



{x e 



X 



1}. Inspired by section 2.3 we have the 



Definition 3.15. Let Hm be a separable Hilbert space of functions defined over Ai and suppose 
Vm is a separable Banach space of functions, also defined over Ai, such that Vm C Hm 
continuously and densely. Let 



A M : [0, T] x n x V M 
B M : [0, T] x Q x V M 



V* 

L 2 (U, H M ) 



be progressively measurable, where U is a fixed separable Hilbert space and Am is a differential 
operator on Ai . Then the equation 



dX{t) = A M {t, X{t)) dt + B M (t, X{t)) dW(t) 
X(0) = x 



(3.2) 



where W(t), t 6 [0, T], is a U -valued cylindrical Q-Wiener process with Q = I is called a 
stochastic partial differential equation on Ai. 



We employ assumption 2.2 on Am and Bm and so the way we define what we mean by a 
solution to (3.2) is 

Definition 3.16. A continuous H M -valued (J 7 ^- adapted process X{t), t 6 [0, T], is called a 
solution of (3.2), if for its dt®¥ -equivalence class X we have X G L a ([0,T] x fi, dt <g)P; Vm) H 
L 2 ([0,T},n,dt ®F;H M ) with a as in H3 andF-a.s 



X(t)=X(0)+ / A M (s,X(s))ds 



B M (s,X(s))dW(s),te[0,T], 

where X is any VM-valued progressively measurable dt ® F-version of X. 

This is completely analogous to definition |2.3| replacing V, H, A and B with Vm, Hm, Am and 



Bm respectively and we immediately have from theorem 2.4 
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Theorem 3.17. Let Am and Bm satisfy assumption 2.2 and suppose X Q £ L 2 (Q, J^o, P; Hm)- 
Then there exists a unique solution to (3.2) in the sense of definition 3.16. Moreover, 



E 



sup ||X 

te[o,T] 



2 



< OO. 



We see that the abstract theory of SPDEs on Ai is a special case of the abstract theory of 
SPDEs, established in chapter 2. We proceed to show that the abstract objects Vm, Hm,Am 
and B_m actually exist, by giving two examples. 



3.3 The stochastic heat equation on a Riemannian manifold 

Here we generalise section [2~3] to Ai, where Ai satisfies assumption |3.14| Let Am := Am, the 
Laplace-Beltrami operator on Ai. Recall from (3.1) that 

A M u = 



)=d m (V\g\g mk dku) 



\9 

in local coordinates, where Einstein summation notation is used. 

We proceed to define the following Lebesgue and Sobolev spaces as given in Hebey 1996 
page 10. 



Definition 3.18. We define the norms 



\u\\lp(M) 



\U\\wi, P (M) 



\u 



M 



p Ma)) 



\ i/p 



HI 



LP{M) + II^ M IIlp(X) 



1 < p < OO 
1/p 



1 < p < OO 



where IIVwl 



Lp(M) 



Vm| ||ip(A^) an d is the covariant derivative of u with (Vu) 



diu in 



local coordinates. 

We define, for 1 < p < oo the spaces 



L P (M) := {u E C°°(M) : \\u\\ LP{M ) < oo} 



W^iM) := {u e C°°{Ai) : \\u\\ w ,, v{M ) < oo} 



\\lp(m) 



Wo' p {M) := {u e C™{Ai)} 



II ■ llivl>P(A<) 



where C^°(Ai) is the space of C°°(Ai) functions with compact support. For p = 2 we use the 
notation of 

H\At) = W^ 2 {M) 
Hl{M) = W^ 2 (Ai). 

The notation o/c" ° means the completion of space C with respect to the D-norm. 



We proceed to briefly discuss Sobolev embeddings for the above spaces. We follow Hebey 



[1996 2000 for the following discussion. 

Recall from when A is an open and bounded subset of M n that Hq (A) ^ -f^ 1 (A) for non-zero 
constant functions are in H l (A) but not in Hq(A). However, when Ai is complete (as in our 
case) we have that (Hebey |1996|, theorem 2.7) 



Wq' p (M) = W l ' p (M) for all p > 1. 
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Thus in our case we have Hq(M) = H 1 (A4). 

Furthermore, the Rellich-Kondrakov theorem for open bounded subsets of M n (Adams 
|2003j) is generalised to the Ai that we are considering via (Hebey 2000 theorem 2.9) 



Theorem 3.19. Let Ai be a Riemannian manifold satisfying assumption 3.14 



(i) For any q £ [0, n) and any p > 1 such that 1/p > 1/q — 1/n the embedding ofW 1 ' q (A4) 
in L p (At) is compact. 

(ii) For any q > n, the embedding of W l,q (A4) in C°(Ai) is compact. 

Remark 3.20. Some comments are needed on theorem \3.19\ 

(i) First of all, the full generality of the theorem has not been stated. For the general statement 
and proof the reader is directed to Hebey 1 2000^ page 37. 



(ii) From part (i) of the theorem, one can choose p = q to see that W 1,9 (Ai) CC L q (Ai) for 
every 1 < q < n. 

(Hi) From part (ii) of the theorem, we see that W 1,q (Ai) CC L q (Ai) for any q > n. Indeed, 



this follows as C°(A4) C L q (Ai) for any q > n. Indeed, by using the arguments of Evans 
199^ one has that 



W 1,9 (M) CC L q (M) for every 1 < q < oo. 



(3.3) 



Finally, we have that the Poincare inequality in an open, bounded subset of R n (Adams 
|2003| ) is generalised to the Ai that we are considering via the following theorem. 



Theorem 3.21. Let Ai be a Riemannian manifold satisfying assumption\3. 14 and let 1 < q < 



oo. Then there exists C p = C p (Ai, q, n) > such that for every u £ W 1,q (Ai) 



u — u\ q du(g) 



M 



< C p [J \Vu\ q dv(g)^ 



where 



Vol(M) j M 



udu(g). 



Proof. Fix 1 < q < oo. Inspired by the analogous proof in the Euclidean case QEvans|[l998l ), 
suppose the above is false. Then we can find a sequence Uk £ W l > q (AA) such that 



\Uk — U k \\Li(M) > k\\ Vu k \\ L q( M ). 



Define 



\ u k — Uk\\Li 

then ||ffc||L9 = 1 and Vk = for every k £ N. Note that llVt^H^ < 1/k and so (vk) is a bounded 
sequence in W l,q (Ai). In light of remark 
v £ L q (M) such that v k] v in L q (M] 



3.20 



there exists a subsequence Vk in W 1,q (Ai) and 
as j —> oo. Thus, by above \\v = 1 and v = 0. 
Since IIVi^Hm < 1/k for every k £ N, we have that v £ W 1,q (Ai) with Vf = a.e. Since Ai 
is connected this implies v is constant. Since v = and v constant this implies that v = and 
so \\v\\iq = which contradicts the above which says that \\v\liq = 1. □ 
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The reader should be aware that the above theorem is only found for 1 < q < n in Hebey 1996 



2000 . Inspecting the proof as given in |Hebey|[l996|[200"0"l , it seems as though this restriction 



of q is due to the method of the proof. 

We see immediately that if u G and J M udv(g) = then 

|M|l 2 (.M) < C p \\Vu\\ L 2( M y 



However, in light of remark 2.7 since we are using the Laplace-Beltrami operator, we will see 
that we do not need to use Poincare, which is advantageous as asking a function to have 
integral may not be what is required in a mathematical model. 

Now take Vm '■ = H l (M) and Hm := L 2 (Ai). Subsequently, we drop the subscript M. for 
the rest of this chapter. Note by definition |3.18| we immediately have the following 

Proposition 3.22. The space C°°{M) is a dense subspace ofV and V C H both continuously 
and densely. Consequently, identifying H* with H, we have the Gelfand triple V C H C V* . 

Up to now, we have only commented on the operator Am- For the operator Bm, let U be 
a separable Hilbert space and let i : U — > H be Hilbert-Schmidt. By proposition 2J3 such i 
exists and so we have now formulated the stochastic heat equation on A4 by 



dX(t) = A M X(t) dt + i dW{t) 
X(0) = X 



(3.4) 



where X e L 2 (Q, Fq, P; H) is given. Note here (fl, J 7 , P) is a complete probability space and 
the cylindrical Wiener process W(t), t e [0, T], is with respect to a normal filtration (F t ) 
analogous to the stochastic heat equation on an open subset of lR ra of section |2.3 



The existence and uniqueness of a solution to (3.4) is covered by the following. 

Theorem 3.23. Let X 6 L 2 (f2, F Q , P; H). Then there exists a unique solution, in the sense 
of definition 3.16, to equation ( 3.4[ ). Moreover, 

E[ sup ||X(t)||y < oo. 
te[o,T] 



Proof. It suffices, by theorem 3.17 , to verify that assumption 2J2 hold for A := and B := i. 
To this end 



1. Since as A is linear HI is satisfied. 

2. To see H2, observe that as B is independent of the solution X we have that for any 
u,v G V, \\B(-,u) - B(-,v)\\ L2 (u,h) = 0. Since C°°(M) is dense in V, for u,v 6 V 
arbitrary, there exists Uk,Vk G C°°(A4) such that Uk — > u and Vk — > v in V as k — > oo. 
Hence, as M. is without boundary 



2(A. M (u -v),u-v) = 2 lim (A M (u k - v k ),u k - v k ) 

k— >oo 



H 



—2 lim ||V(ufc 

k— >oo 



-2\\V{u-v)\\ 2 H < 



Thus H2 is satisfied with c = 0. 
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3. To see H3, using the same argument as above for v G V one has 

2(A M v,v) = -2||Vv||^ = 2\\v\\ 2 H - 2\\v\\ 2 v 

since \\v\\ v = \\v\\%+ ||Vw||#. Recall that as % is Hilbert-Schmidt there exists C4 > such 
that ||i||L 2 (£/,ff) < c 4 an d so 

2(A M v,v) + \\i\\l 2{u>H ) < 2\\v\\% - 2\\v\\v + 4 

Noting that cl is (J r t )-adapted and is in L 1 ([0, T] x Q, dt <8> P) we see that H3 is satisfied 
with a = Ci = c 2 = 2 and f(t) = cf. 

4. Finally, for H4 let u, v G C°°(.M). Then as .M is without boundary 

|(Axix,u)| = |(Axw,u)jt| = |(Vu, Vu)i?| < ||Vu||jf||Vu||jj < ||u||y||u \\v- 

This implies that ||A^u||y* < \\u\\y for all u G V by a density argument and so H4 is 
satisfied with c 3 = 1 and g(t) = 0. 



We now apply theorem 3.17 to see that (3.4) has a unique solution. □ 



3.4 A nonlinear stochastic partial differential equation on a Rieman- 
nian manifold 

Until now, we have only considered linear SPDEs. In this final section, we will look at a specific 
nonlinear SPDE. We replace the Laplace-Beltrami operator in the stochastic heat equation with 



the p-Laplace-Beltrami operator where p > 2, which generalises example 4.1.9 of Prevot and 



Rockner 2007 to our manifold Ai. 



Let M. be a Riemannian manifold satisfying assumption |3.14| Define 

V := {u G W l ' p (M) : I udu(g)=0} 

Jm 

H := {u G L 2 (M) : I udu{g) = 0} 
Jm 

where p > 2 and equip V and H with the W l,p (A4) and L 2 (Ai) norms respectively. We see 
that C°°(M) n V is dense in V and V C H continuously and densely. Hence V C H C V* is a 
Gelfand triple. 

Define A : V — y V* by 

Am := div^dVwl^ 2 Vw), 
by which we mean for given u G V, 

(Au,v) '■— — I |Vm| p ~ 2 (Vu,Vv) g di>(g) for every v G V, 
Jm 

where |Vw| p := (g iJ d Xi ud x u) p ^ 2 and (Vw, Vt> ) g := g l W Xi ud x .v. 
For li.oey 



|(Au,v)|</ |V«r _i |Vu| di/(p) < / \Vu\ p dv{g)\ [ \Vv\ p dv{g) 

Jm \Jm J \Jm 

11 iip-in n 
< \\u\\y 
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which implies that 

||Au||v* < for every u G V. (3.5) 

This shows that Au is a well defined element of V* and is bounded as a map from V to V*. 

For the noise term, as before fix U a separable Hilbert space and let W be a [/-valued 
cylindrical Q-Wiener process with Q = I. Let i : U — > H be Hilbert-Schmidt, which by 
proposition 2J3 always exists. 

We now have 

Theorem 3.24. Let Xq G L 2 (f2, J-" 0; ^0- Then there exists a unique solution to 

dX(t) = dw M (\VX(t)\ p ~ 2 VX(t))dt + idW(t) (p > 2) 
X(0) = X , 



in the sense of definition \3.16[ Further, 

E[ sup ||X(t)||y < oo 

«e[o,T] 

Proof. As before, it suffices to check that A and B := i satisfy the hypotheses of HI to H4 of 



assumption 2.2 



1. To check HI it suffices to show that for u, v, x 6 V and A G K with |A| < 1 

lim / (\V(u + Xv)\ p ~ 2 (V(u + Xv),Vx) g -\Vu\ p - 2 (Vu,Vx) g ) du(g) = 0. (3.6) 

J M 

Clearly the integrand converges to zero as A — > 0, so we need only find an L l (M) bounding 
function (independent of A) to use Lebesgue's dominated convergence theorem. 

To this end, since |A| < 1, using Cauchy-Schwartz and the fact that x h- > x q is convex for 
q > 1 one immediately has 

| \V(u + Xv)\ p ~ 2 (V(u + Xv), Vx) g \ < 2 p - 2 (\Vu\ p - 1 + iVvf' 1 ) \Vx\ 

and so the integrand is bounded above by 

(2 p - 2 + 1) IV^ 1 |Vx| + 2 P ~ 2 l^vf- 1 \Vx\ 

which is clearly in L l (Ai) and so applying Lebesgue's dominated convergence theorem, 
we see that (3.6) follows. 

2. For H2, since B is independent of the solution, for u,v G V it follows that \\B(-,u) — 
B{- ) v)\\l 2 (u,h) = 0- Further, using the Cauchy-Schwartz inequality one has 

-(Au- Av,u - v) = / (|Vw| p ~ 2 Vw - \Vv\ p ~ 2 Vv, Vu - Vv} g du(g) 

Jm 

I \Vu\ p - \Vuf~ 1 \Vv\ - iVvf' 1 |V«| + \Vv\ p dv(g) 
Jm 



> 



[ {iVur 1 - iVi^r" 1 ) (|Vu| - |Vv|) du(g) 
Jm 



> 

where the last inequality holds since s i— > s q is increasing for s > and q > 1. Thus H2 
holds with c = 0. 
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3. To see H3, using the Poincare inequality (theorem 3.21) and the definition of V there 



exists C p > such that 



/ \Vu\ p du(g) > C p 1 / \u\ p dv(g) for every u G V 
Jm Jm 



IM JM 
and so for all v G V 

(Av,v) = -\\Vv\\ p LP{M) < -C- 1 ||«||^+ ||Vv||^, (M) - \\Vv\\ p LP{M) 
< max(-l, -Cp^Wvfy + \\Vv\\ p LP{M) 

Thus 



which implies 



-\\Vv\\ p LP{M) < -min(l,C p - 1 )|| U ||^ + \\Vv\\ p LP{M) 



(Av,v) = -\\Vv\\ p < ^r-t-t-M* 



\LP(M) — 2 n^iiv- 

Since i : U — > H is Hilbert-Schmidt, there exists C4 G (0, 00) such that ||«||l 2 (c/,h) < C4, 
thus 

2(Av,v) + \\i\\l 2{UjH) < -wmfrC^WvWlr + 4 

which shows that H3 is satisfied with a = p > 2, c\ = 0, C2 = min(l,C ? 7 1 ) > and 
f{t) = c\. 



4. Finally, H4 follows from (3.5) 



Thus applying theorem 3.17| completes the proof. □ 



4 Stochastic partial differential equations on moving sur- 
faces 

4.1 The stochastic heat equation on a general moving surface 

In order to build up intuition as to what a SPDE on a moving surface should look like, we first 
consider the deterministic case. 

4.1.1 The deterministic case 

Let A4(t) be a hypersurface for each time t G [0,T] where T G (0, 00) is fixed. We need some 
notion of what it means to have such an object. Unless otherwise stated, the definitions and 



proofs are found in Deckelnick et al. 2005 



Definition 4.1. Let k G N. A subset F C IR n+1 is called a C k -hypersurface if for each point 
xq G T there exists an open set U C IR n+1 containing xq and a function <p G C h (U) such that 

U n r = {x G U I (j){x) = 0} and V0(x) ^ for every x E U (IT. 

This allows us to define what it means for a function on T to be differentiable. 
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Definition 4.2. Let F C IR n+1 be a C 1 -hyper surf ace, x G V . A function f : F — > R is called 
differentiable at x if f o X is differentiable at X _1 (x) for each parameterisation X : — > R n+1 
ofF with x G X(Q). 

The following lemma shows us how to interpret the above definition in terms of functions 
defined on the ambient space. 

Lemma 4.3. Let F C R n+1 be a C 1 -hypersurface with x G F. A function f : F — > R is 
differentiable at x if and only if there exists an open neighbourhood U in IR n+1 and a function 
f : U — > R which is differentiable at x and satisfies f\rnu — f ■ 

With the notion of differentiable functions on F we now define the tangential gradient, which 
is the form of the differential operator we will be considering. 

Definition 4.4. Let F C R ra+1 be a C 1 -hypersurface, x G F and f : F — > R differentiable at x. 
We define the tangential gradient of f at x by 



V r /(x) := V/>) - (v/(z) ■ v{xj) u(x) 



Here f is as in lemma 4-3, V denotes the usual gradient in M. n+1 and u(x) is a unit normal at 
x. 

This leads to the definition of the Laplace-Beltrami operator on F(t), 

Ar(t)f '■— Vr(t) • Vr(j)/. 

In the following let X G C 2 (M(0) x [0,T]; IR n+1 ) be a local parameterisation of M(t), where 
we assume that M.{t) is compact, connected, without boundary and oriented for every t G [0, T}. 
We assume that points on Ai(t) evolve according to X t (x, t) = V u , x G .M(O) where V u is the 
velocity in the normal direction and that X(-,t) : Ai(0) — > M.{t) is a diffeomorphism. We 
define the Sobolev spaces H l (Ai(0)) and L 2 (Ai(0)) with respective norms analogously as given 
in definition 13.181 

Before stating the conservation law and deriving the PDE, we need to define a time derivative 
that takes into account the evolution of the surface, generalise integration by parts and give 
the so-called transport theorem. 

Definition 4.5. Suppose F(t) is evolving with normal velocity v v . Define the material velocity 
field v := v v + v T where v T is the tangential velocity field. The material derivative of a scalar 
function f = f(x,t) defined on Qt '■= U t ^[o,T]F{t) x {t} is given as 

df 

We now give a generalisation of integration by parts for a hypersurface F, the proof of which 



is found in Gilbarg and Trudinger 2001 



Theorem 4.6. Let F be a compact C 2 -hypersurface with boundary and f G H /1,1 (r; R n+1 ). 
Then 



V r • fdH n = f- HvdH n + / / • v 8v dH 
r Jt Jar 

where H = Vr • v is the mean curvature and vqy is the co-normal. 



n-l 
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This leads us nicely onto the following lemma which is referred to as the transport theorem, 
whose proof is given in Dziuk and Elliott 2007 . 



Lemma 4.7. Let C(t) be an evolving surface portion ofT(t) with normal velocity v v . Let v T be 
a tangential velocity field on C(t). Let the boundary dC(t) evolve with the velocity v = v v + v T . 
Assume that f is a function such that all the following quantities exist. Then 

d 
d7 



/= / d'f + fVr-v. 

C(t) JC(t) 



We now have all the necessary theory to formulate an advection-diffusion equation from the 
following conservation law. 

Let u be the density of a scalar quantity on T(t) and suppose there is a surface flux q. 
Consider an arbitrary portion C(t) of T(t), which is the image of a portion C(0) of T(0), evolving 
with the prescribed velocity v v . The law is that, for every C(t), 

d 
df 



u = - q-Vdv- (4.1) 

C(t) JdC(t) 



Observing that components of q normal to C(t) do not contribute to the flux, we may assume 
that q is a tangent vector. With this assumption, theorem 4J3 lemma 4/7 and assuming 
q = uv T — Vr(t)M one has the PDE 

d'u + uVr(t) • v — Ar(t)it = 0. 

We now take T(t) = Ai(t) and assume for simplicity that v T = 0. In this case we have that 
v = Vv and so VmU) ' (V 11 ) = VH, where H is the mean curvature of Ai(t). We now arrive at 
the following model PDE on Ai(t) 



(4.2) 



d'u + uVH - A M( t)U = 

u(x, 0) = u x e M(0). 
For x G Ai(0) define w(x,t) := u(X(x, t), t). Then w is defined on Ai(0) and 

dw Ou 

— (x, t) = ~^{X{x, t),t) + (Wu)(X(x, t),t)- X t (x, t) 

= ^(X(x,t),t) + (Vu)(X(x,t),t) ■ V u (X(x,t)) (4 ' 3) 
=: O'u. 



Further by Deckelnick et al. 2005 , letting y = X(x,t) one has 

AMt " u{y) = vikwi k i 9 "^ (xA (44) 

where gij(x,t) := X x .(x,t) ■ X x .(x,t), g^(x,t) is the (i,j) th element of the inverse of g(x,t) : = 
(gij(x, t))ij and \g(x, t)\ := det g(x, t). We employ the Einstein summation notation and assume 
that there exists k\ > such that \VH\ < k\ for any (x,t) G M(t) x [0, T]. Note here that 
Xj are not the local coordinates of x, but are the local coordinates of a parameterisation that 
gives x. 
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Putting all this together, we see that w solves 
^(x,t) +w(x,t)V(X(x,t))H(X(x,t)) - J_ A (V(M)v^M^') (x,t) = 

w(x, 0) = M 

(4.5) 

which we solve on M.(0). On solving, we set u(y,t) := w(X 1 (y,t),t). 

We will drop the V(X(x, t))H(X(x, t)) notation and simply write VH in the following. We 
see that we have reduced the PDE on a moving surface to a PDE on a fixed surface, M(0). 
This will allow us to define the stochastic analogue, but importantly we must define what noise 
we are considering. 

Remark 4.8. Since det^DX^ 1 (■,{)) is continuous, bounded and bounded away from for every 
t G [0,T] there exists ai,b\ > such that a\ < det(DX _1 (x, £)) < b\ for every (x,t) G 
M.{t) x [0, T]. By the smoothness of the parameterisation and the compactness of M.(fS) x [0,T] 
there exists a-i, &2 > stzc/i 



a 2 



< y/\g(x,t)\ < b 2 for every (x,t) G M(0) x [0, T]. 



Furthermore, since (gij(x,t))ij is positive definite and symmetric for every (x,t) G Ai(0) x 
[0,T] it follows that (g^ (x,t))ij is also positive definite and symmetric and so since {g^)%j 
contains functions which are continuous and M.(0) x [0, T] is compact there exists 03,63 > 
such that 

~ 01) (jV ~ 

a 3 \Vv\ 2 < g %3 {x,t) < h\Vv\ 2 for every v G H^M^O)), (x, t) G M(0) x [0,T] 

where |Vt>| 2 is notation for Y17=i(^i v ) 2 an ^ ^ ^ s no ^ ^ e 9 ra dient on the ambient space. 
Finally, by the compactness of M.(0) x [0,T] there exists b > such that 

\g ij (x,t)d x .ud Xi v\ < b\Vu\ \Vv\ for every u,v G H\M{0)), (x,t) G M(0) x [0,T]. 
4.1.2 The stochastic case 

Let H := L 2 (.M(0), dv(^o); R) and fix f/ a separable Hilbert space. Here di/(g ) is the Rieman- 
nian volume element, not to be confused with the unit normal, v. Define 

H gt := L\M(0),^/\gi^t)\du(g o y,R). 



By remark 4.8 it is easy to see that H and H gt coincide as sets since 

02|MIh < h\\ 2 Hgt < b2\\u\\ 2 H . (4.6) 

Let If be a U— valued cylindrical Q— Wiener process with Q = I. Let i : U — > H be Hilbert- 
Schmidt, noting that proposition |2.5| ensures that such i always exists. Define 



G t : H ^ L 2 (M(t),du(g t );R 

by 

(G t f)(x) := f(X-\x,t)) 
noting that the following shows that this map is well defined. 
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Lemma 4.9. Suppose that v(-) G H. Then v(X 1 (-,t),t) G L 2 (M(t),dis(g t );M) for every 

te[o,T]. 

Proof. We see that there exists c > such that 



c> \v\ du(g ) 

Im{o) 



and so using remark 4.8 we have 



c > I \v\ 2 du(g ) > ai a 2 I \v(X~\y,t),t)\ z dy > ^ / \v(X-\; t)) | 2 dv{g t ) 

JM(0) JR" °2 JM{t) 

which completes the proof. □ 
We define the noise on Ai(t) by idW x 1 ('' t \t) which is defined as 

idW x ' 1{ -' t \t) := G t idW(t) (4.7) 

and the above shows that the noise is L 2 (A4(t), dh>(g t ); M.) valued. 
We now define the stochastic analogue of (4.2) as 



d'u = (A M{t) - VH)udt + idW x 1{ -' t} (t) 
u(0) = u 



(4.8) 



which we interpret as solving the following SPDE on M.(0) (in the sense of definition 3.16) 
with Gelfand triple H\M(0)) C L 2 (M(0), y/\g(-,t)\ dv(g ); R) C (H^MiO)))* 



d^= f ^_^ ^(^(x,t)y/\^^\p^)-VH)wdt + idW(t) 
\y/\g(x,t)\dxj V dxj J 



Note here that the operator is in local coordinates here and that i dW(t) is if- valued noise, but 
by (4.6) we see that it is if 9t -valuecj^] This gives 



Definition 4.10. Suppose that we can solve equation (4.9). Call the solution w. Then we 
define the solution to equation (4.8), u by 

u(t,u)(y) :=w{t,u){X-\y,t)) 



where lo G fi, t G [0,T] and y G M.(t). Here we adopt the notion of solution to (4.9) in the 
sense of definition 3.16] 

From the definition of L 2 (M(0)) it follows that H 1 (M(0)) C L 2 (M(0)) continuously 
and densely and so by the equivalence of the H and H gt norms we have that if 1 (A^(0)) C 
L 2 (A4(0), a/|(?(-, t)\ dv(g ); M) continuously and densely and so indeed 

H\M(o)) c L 2 (M(o),^A^W\Mgo);^) c (H^Mio)))* 

is a Gelfand triple. 

For brevity, we let V = if^A^O)) and H g = L 2 (M(0), y/\g(;t)\ du(g ) ; R) (so we drop the 
subscript t). The following shows that we can solve (4.9). 



2 Strictly speaking, we should replace i by tp t i where tp t : H — > H gt is given by <p t f(x) — f(x)/(\g(x, t)\) 1 / A . 
Then ||<^t||o P = 1 and we may consider i : U —> H gt as Hilbert-Schmidt. 
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Proposition 4.11. Suppose uq G L 2 (Q, J-'o, P; H). Then there exists a unique solution of (4.9), 
in the sense of definition \3.16 . Moreover, the solution w satisfies 

E[ sup ] < oo. 

te[o,T] 9 



Proof. By theorem 3.17 it suffices to show that 

1 d 



A :-- 



B 



V\g{x,t)\dxi 



^(x,t)^\gjx~t)\ 



_d_ 

dx,, 



VH 



satisfy HI to H4 of assumption |2.2| To this end 

1. Clearly as A is linear, HI is satisfied. 

2. For H2, we use the pairing defined by (z,v) 9 = (z,v)h s for every z G H g , v G V 
defined in the obvious way. Let u, v G V and since B is independent of the solution we 



have that \\B(-,u) — B(-,v)\\L 2 (u,H g ) = 0. By the arguments of the proof of theorem 3.23 



we see that integration by parts is valid for elements of V and we see that we can identify 
the pairing of V and V* with the inner product on H g . Hence 



{A(u — v), u — v) g 



M(o) 



d Xi {g %3 {x,t)y/\g(x,t)\d Xj {u -v))(u-v) du(g ) 
y/\g(x,t)\VH(u-v) 2 dv(g ) 



< 



M(0) 

a 2 



\u — v\ 



where the inequality follows from the positive definiteness of and the equivalence of 
the H and H a norms. Thus H2 is satisfied with c = > q. 

9 a 2 



3. For H3, let v G V and fix t G [0, T\. Then using remark 4.8 and noting that v is time- 
independent one has that 



(Av,vY 



g 13 (x, t) yJ\g{x,t)\d Xj vd x .v dv{g ) 

M(0) 

y/\g(x,t)\VHv 2 du(g ) 



>M(0) 

< -a 2 a 3 



< 



a 2 a 3 



M{0) 

f 

M(0) 



\Vv\ 2 dv(g ) + b 2 k 1 \\v\\ 2 H 
g l: >(x } 0)d Xt vd x „vdis(g ) + b 2 ki\\v\\ 2 H . 



I I Z 

Now identifying that g %1 {x,{))d x .vd Xj v = \^m(o) v \ we have 



(Av,v) 9 < 



a 2 a 3 



Vm(o)v\\ 2 h + b 2 k 1 \\v\\ 2 H 



< hh + 



a 2 a 3 



h ) 



a 2 a 3 



\H 



IV) 
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where the last inequality follows by the definition of the V^-norm. Now using the equiva- 



lence of the H and H g we have 



2(Av,v) 9 + \\i\\ 2 LmHg) < Cl \ 



a , 2 



C2\\V\\ V + C 4 



where a = 2, c x = 2(^ + g) > 0, c 2 = ^ > and c 4 > \\i\\ L2 (u,H a ) with c 4 existing 
and finite as % is Hilbert-Schmidtj^J which shows H3. 



4. Finally for H4, let u,v G C°°(^W(0)) and again using remark 4.8 and noting that u and 
u are time independent one has that 



\(Au,v) 9 \ < 



g tJ d x ud x .vdv(g 



+ b 2 ki\\u\\ H \\v\\ H 

1/2 



'M(O) 

<b 2 b( [ \Vu\ 2 du{g 

\JM(0) 

b 2 b 

< I|V^(o)M||h||V>i(o)V||^ + &2&l|MI# 

&2 



M(0) 



1/2 



V« I du(g ) + 6 2 A;i||u||h||v||h 



< 2 max 



hb 

a 2 ' 



b 2 k\ \\u\\y\\v\ 



v 



which implies that ||Au||v« < C3||w||v which, by a density argument, gives H4 with C3 
2 max (g,M:i) > and g(t) = 0. 



□ 



The following gives a regularity estimate for the solution u of (4.8). 
Proposition 4.12. Suppose u$ G L 2 (Q, J-q, P; H g ). Then the solution u of (4.8) satisfies 



E[ sup ||w(t)||i2(^( t ))] < 00. 
te[o,T] 



Proof. From proposition 4.11 and the equivalence of the H and H q norms one has that 



a 2 E[ sup |K*)||y <E[ sup |K*)|||-] < 00, 
te[o,T] ' te[o,T] 



where w is the solution to (4.9). However, by lemma 4.9 we have that 



W(t)\\ 2 L 2 {Mm = \\w(X 1 (- 1 t))\\ 2 L 2 {M{t)) < —\\w u ' ,i: ' 



a-2 



H 



and so taking the supremum over all t G [0, T] and then taking expectations yields the result. □ 



Remark 4.13. There still remains the question of uniqueness of the solution u to (4.8). There 



was a choice of diffeomorphism to take and we always ensured that the parameterisation and 



the diffeomorphism were compatible, in the sense that (4.3) holds. Thus, we only speak about 
uniqueness up to parameterisation. 



3 When i is considered as in the footnote 2. 



23 



4.2 The stochastic heat equation on a sphere evolving under mean 
curvature flow 

We now give a specific choice of -M(i), namely the 5"™ _1 sphere evolving under so called 'mean 
curvature flow'. 

Definition 4.14. Let T be a C 1 hypersurface with normal vector v . We define the mean 
curvature at x G V as 

H(x) := V r • v. 
This naturally leads us onto the following 

Definition 4.15. Let (r(t)) te [ 0T ] be a family of hyper surf aces. We say that T(t) evolves ac- 
cording to mean curvature flow (mcf) if the normal velocity component V satisfies 

V = -H. 



For our case, as given in Deckelnick et al. 2005 , one defines the level set function by 



(x,t) = ||x|| — R(t), which describes a sphere of radius R(t). Indeed, by Deckelnick et al. 



2005 , one has 

n 

H = V • = -, 

where V is the gradient in the ambient space. Further, V = <p t = —R- Hence solving V = —H 
yields R(t) = y/1 — 2nt for t G [0, ^-), noting that the initial radius is 1. We observe that at 
t = l/2n the sphere shrinks to a point and so for the remainder for this section we will fix 
T < l/2n. 

From this, we see that we will consider 

S(t) := {x G R n : ||x|| = y/l - 2nt} t G [0,T]. 

Observe that <S'(0) = S^ 1 . Indeed with this representation of S(t) we have the following natural 
parameterisation and diffeomorphism 

X(-,t) : S(0) — > S(t) x ^ X(x, t) := xy/l - 2nt. 

We then see that 

gij(x,t) = (1 - 2nt)gij(x,0) 

and so as gij(x,0) is diagonal 



We now use the PDE (4.2), which yields 

,2 



71 

d ' U -T^2n-t U - As ^ U = ° (4.10) 
u(x, 0) = Uq x G 5(0). 



Letting w(x,t) = u(X(x,t),t) as done in section[4X yields the PDE on S(Q) as 

dw . , n 2 1 d ( ,•„•. . /-— -rdw\ , 

— (x,t) - w(x,t) — g l Hx,t)J \g(x,t) — (x,t) = 

dt K ' 1 >l-2nt y/\g{x,t)\dXi \ y 1 >V ^ dxJ K ' (4.11) 

w(x, 0) = Mo 
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which we solve on 5(0). On solving, we set u(y,t) := W7(X _1 (?/, t), t). 

By the isotropic nature of the evolution of S(t), we can work without the weighted L 2 (S(0)) 
space. 

For the noise, analogous to section 4.1, we define H = L 2 (S(0)) and let U be a fixed 
separable Hilbert space. Let i : U — > H be Hilbert-Schmidt, which by proposition 2.5| always 
exists. Let W be a ^/-valued cylindrical Q- Wiener process with Q = I. We define the noise on 
S(t) by 

zdW^ 1 ^ 



which is given by (4.7). We define the stochastic analogue of (4.10) as 



d'u = (A 5W + u dt + i dW*- 1 ^ 

u(0) = u 



(4.12) 



which we interpret as the following SPDE on S(0) (in the sense of definition 3.16) with Gelfand 
triple H\S(0)) C L 2 (S(0)) C (^(5(0)))* 



dw 







\f\9{x,t)\dxj 



f\x,ty\g{x,t)\ 



dxij 



+ 



n 2 w 



2nt 



dt + idW{t) 



(4.13) 



w(0) = Uq 



We define the solution to (4.12) analogously as in definition 4.10, namely 

u{t,u){y) :=w(t,u)(X- 1 (y,t)). 



The following shows that we can solve (4.13). 



Proposition 4.16. There exists a unique solution to (4.13) in the sense of definition 3.16, 
Moreover, 

E[ sup imt)||#] < oo 
*e[o,T] 



Proof. By theorem 3.17 it suffices to show that 

1 d 



A :-- 



y/\g{x,t)\dXj 



i ^x,t) v ^gJx~W\ 



d_ 

dx,; 



+ 



n 



1 - 2nt 



B := i 



satisfy HI to H4 of assumption 2.2 In the following, we use the arguments of the proof of 
theorem |3.23| to justify the integration by parts and identifying the pairing between V and V* 
with the inner product on H. 

1. Clearly, as A is linear, HI is immediately satisfied. 

2. For H2, let u, v G V. Noting that \\B(-,u) — B(-,v)\\l 2 (u,h) = and using the isotropic 
evolution of S(t) one has 



(A(u — v),(u — v)) 



1 



1 - 2nt 

■2 



5(0) 



dxj 



dx.-,. 



+ 



rr 



1 - 2nt 1 
1 



\u — v\ 



H 



\-2nV 

,2 



7 s(o){u - v)\\ 2 H + - ™ 2nt \\u - v\\ 2 H 



< 



n 



2nT 



\u — v\ 



H ■ 
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Hence H2 is satisfied with c : = 



2/r 



\-2nT 



> 0. 



3. To see H3, let v G V and then by item 2 above, 

1 



(Av,v) 



n 



Vs(o)u||h + ~ t — ii' ii.//- 



1 - 2nt" " " n 1 - 2nt 
However, since T < l/2n, we have that for every < t < T 

1 1 

1 < < 



1 - 2nt ~ 1 - 2nT 



so 



(A-Ujv) < -||V S ( )f HI + - — 



1 + 



by definition of the norm on V. Hence 



n 



1 - 2nT 



1 - 2nT 1 



\H 



IV) 



2(Au,v) + ||i| 



L 2 (U,H) 



— c l\\ v \ 



H 



c 2 Hy + c : 



where a = 2, ci = 2 ^1 + jr^^j > 0, C2 = 2 and C4 > ||i||.L 2 which exists as i is 
Hilbert- Schmidt. This shows that H3 is satisfied. 

4. Finally, for H4 let m,d6 C°°(S(0)). Then by the isotropic evolution of S(t) we have 

1 



\(Au,v)\ < 



1 - 2nt 
1 

1 - 2nt 
1 



5(0) 



^'M)|^d% ) I 



(Vs(o)tt, V s( o)w) 9 di/(^o) 



1 - 2nT 

2 



rr 



H IIw||h-|M|w- 

1 -2nT H 



'X(O) 

2n 2 



< 



1 - 2nT 



U VP V 



which shows that ||v4u||y* < rz^r lKllv which, by a density argument, shows that H4 
holds with c 3 = ^ nT > and g(t) =0. 



□ 



Analogous to proposition |4.12| we immediately see that 

E[ sup ||tt(*)||£a(5(t))] < 00 
te[o,T\ 



where u is the solution to (4.12). 
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4.3 A nonlinear stochastic heat equation on a general moving sur- 
face 

So far in this chapter, we have only considered linear SPDE. Since some mathematical models 
need nonlinear terms to be more realistic, we present an example of a nonlinear SPDE. 
Recall section 4.1, but instead of ( 4.2[ ) we consider 



d*u + uVH + f(u) - A M{t) u = 

u(x, 0) = u {x) x G M(0). 



(4.14) 



We will specify how / should behave shortly. As before, we assume that there exists k\ > 
such that \VH\ < k x for every (y,t) G M(t) x [0,T]. 

Using the method of section 4.1 by defining w(x,t) := u(X(x,t),t) where x G M(0), one 
immediately arrives at the following PDE on Ai(0) 



l"(x,t)s/\g(x,t)\ 



dw 
dxj 







(4.15) 



w(x, 0) = uq(x) x G A^(0). 



We define the stochastic analogue of (4.14) as 



d'u = (A M(t yu - f(u) - uVH) dt + % dW x ~ 1( --*\t) 
u(0) = u 



(4.16) 



(where idW x (''*)(£) is given by (4.7)) which we interpret as solving the following SPDE on 
M(0) 



dw(t) -- 
w(0) = Uq 



d 



i \x,t)j\g(^)\ 



V\g{x~t)\dx 
with Gelfand triple V C H g C V* where 



dw 
dxj 



-f(w)-wVH dt + idW(t) 



(4.17) 



V 



H. 



and we define 



9 

H 



H\M(0)), 



L\M(0),V\g(;t)\dv(jg 
L 2 (M(0),du(g );R) 



as in section 4.1 The definition of the solution to (4.16) is as given in definition 4.10 Note 
here that idW(t) is L 2 (A4(0), du(g ); E) valued but by (4.6) we see that it is H g valued. 
We now employ the following assumptions on /. 



Assumption 4.17. Consider (4.17). We assume that f : V — > V satisfies 
(i) f is monotone increasing on H g . That is; for any u,v G V 

(f(u)-f(v),u-v) Bg >0 

(ii) f is Lipschitz on H g and f(0) = 0. That is; there exists cn p > such that 

\\f(u) - f(y)\\ H g < c hp \\u-v\\ Hg for every it, v G V. 
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Remark 4.18. The assumptions are very natural and are the sort of assumptions one finds in 
deterministic PDE theory. 

Note that (ii) above implies that f is continuous on H g . 



One can see that (4.17) is simply (4.9) but with an additional nonlinear operator / : V — » V . 
In light of this observation, the following lemma will save needless repetition. 

Lemma 4.19. Let V C H g C V* be as above. Suppose 



A :-- 



B := i 



1 



d 



j\g{x,t)\dxi 



*(x,t)y/\faj)\ 



_d_ 

dx. 



VH 



satisfy HI to H4 of assumption 2.2, with a = 2 in H3 and \\i\\L 2 (u,Hg) < c 4- Then 

A:=A-f 

B = i 



where f : V — > V and satisfies assumption J^.ll also satisfy HI to H4 of assumption 2.2 . 

Proof. 1. For HI noting remark 4.18 we have that / is continuous on H g . For u,v,x G V 
one has 

(A(u + Xv), x) 9 = (A(u + Xv), x) 9 - (f(u + Xv), x) 9 

and we have that A i— > (A(u + Xv),x) 9 is continuous by assumption. We now see that 
A i — y {f{u + Xv),x) 9 is continuous as (f(u + Xv),x) 9 is a composition of continuous 
operators and so continuous. 

2. For H2, let u,v G V. Then 

{Au — Av, u — v) 9 = (Au — Av, u — v) 9 — (f(u) — f(v),u — v) Hg 
< (Au - Av,u- v) 9 



by monotonicity of / on H g and noting that \/\g(-,t)\ > 0. Hence as B(-, u) 
have 

2(Au - Av,u - v) 9 < 2(Au — Av, u — v) 9 < c\\u - v\ 
by assumption on A. 



B(-, v) we 



l 2 



3. For H3, let v G V. Then 



(Av,v) 9 < (Av,v) 9 + \(f(v),v) Hg \ 

<(Av,vy + h\f(v)\\ Hg \\v\\ Hg 



(4.18) 
(4.19) 



where the last inequality follows from the equivalence of the H and H g norms, (4.6). By 
the Lipschitz property of / on H g and the assumption that /(0) =0 we have 

= ll/(«) - /(0) + /(0)|U S < cnpllnll^ 

and so 



(Av,v) 9 < (Av,v) 9 + 



cii P b 2 



thus under the assumption of the lemma 



2(Av,v) 9 +\\i\\l 2(UiHg) < ci + 



cnpb 2 



C2\\v\\v + c\. 
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4. Finally for H4, note that since there exists c 3 > such that 

||Au||y* < c 3 ||w||v for every u EV 
we have that for u, v G V arbitrary 

\(Au, v) 9 \ < cj\\u\\v\\v\\v- 
Bearing this in mind, one computes 

\(Au,v) 9 \ < \(Au,v) 9 \ + \(f(u),v) Hg \ 

^ C\i P b 2 

a 2 

, ( . c Up b 2 \ 



which shows that 



and so H4 is satisfied. 



Clip&2 



\\Au\\ v * <[c 3 + ^ \\u\\ v 



Immediately we have the following result. 



□ 



Theorem 4.20. Let uo G L 2 (Q, Tq, P; H). Then there exists a unique solution (in the sense of 
definition 3.16) to (4.17). Further, 

E[ sup IKOIIL] < 00 

te[o,T] 



Proof. Recall proposition |4.11| where we saw that 

1 d 



A :-. 



y/\g(x,t)\dxi 



^(x,tW\g(x,t)\ 



d 

dXj 



- VH 



B := i 



satisfy HI to H4 of assumption |2 . 2| with a = 2 in H3. By theorem |3.17| we are required to show 
that 

A:=A-f 
B = i 



satisfy HI to H4 of assumption 2.2 However, we now apply lemma 4.19 to see this and the 
proof is complete. □ 



Remark 4.21. (i) We see that we've shown there exists a solution to (4.16), which is unique 
up to parameterisation and diffeomorphism X . 



(ii) Analogously to proposition Jf.,12 we immediately see that the solution u to (4.16) satisfies 

E[ sup \\u{t)\\ L 2 iM{t)) ] < 00. 
te[o,T] 
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5 Stochastic partial differential equations on general evolv- 
ing manifolds 

5.1 Discussion 

We proceed to give a different (but under some conditions) equivalent way of thinking of an 
evolving manifold. The idea is to think of one fixed topological manifold, Ai, equipped with a 
one-parameter family of metrics {gij{',t))te[Q,T\ applied to the manifold. 

This approach of thinking of evolving manifolds is far from new. It is the standard view 
when one considers Ricci flow of manifolds, for example (Topping 2006| ). 

The idea is now to put a PDE on M. with the metric (gij(-, t))te[o,T\ and define the stochastic 
analogue of this. 

In the following we discuss how, under certain regularity assumptions on the metric, this is 
equivalent to the PDE considered in chapter |4} We then proceed to define what PDE we will be 
considering on M. and formulate the stochastic analogue, proving an existence and uniqueness 
result. 

We will always consider M. to be a compact, connected, oriented and closed topological 
manifold. We also assume that no topology changes occur over [0, T\. 
For the discussion, suppose we are given the metric 

9ij(x,t) = f(t)gij(x,0) x E M, 

where gij(x, 0) is sufficiently nice and / G C 1 ([0, T}; (0, oo)). In order to compare equations in 
this case, we need to find a parameterisation X such that 

9ij{x,t) =X Xi (x,t) -X Xj (x,t), (5.1) 

subject to X t ■ v — v u where v„ is the velocity in the normal direction. 

If the metric is initially diagonal, then it is diagonal for all times and so solving (5.1) is 
equivalent to solving the eikonal equation on Ai. Further, in a special case when gu(x, 0) = 1, 



the existence of solutions are discussed in Kupeli 1995 



Supposing that we can solve (5.1), consider the following PDE on M. 



du 

m- AMU = (5.2) 



u(x, 0) = uq(x) x G M, 
where M. is equipped with the metric {gij{-, £))te[o,T] and so 



in local coordinates. As in section 4.1, let w(X(x,t),t) = u(x,t), where X(-,t) : M. — > Ai(t) 
is the parameterisation which gives rise to the given metric. Then noting that X t (x,t) =: 
v(X(x,t),t) and X t - v = v v one has 

^ = ^(w(X(x,t),t)) = ^(X(x,t),t) + (Vw)(X(x,t),t) ■ X t = d'w 



and 

A . 11 f \ = — i ( , , ! 

\g(x,t)\dxi\ ' ' dxj 



1 d ( du 
&M(t)w{y,t)= — — (g v (x,t)^/\g(x,t)\— ) = A M u. 
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This shows that 

du 

— A M u = onAl 

at 

implies 

d'w — A M ^w = on Ai(t) 



which recalling (4.2), is almost the heat equation on Ai(t). 
Now consider 

du TT . 

- + uH-A M u = (53) 

u(x, 0) = u (x) x G Ai, 

where H is the mean curvature of Ai under the given metric (g%j(-, t))te[o,T\- Then, immediately 
the above discussion shows that we have the following PDE on Ai(t) 

d'w + wH — A M ^w = 

which is the heat equation on Ai(t) if Ai(t) is a hypersurface with evolution completely in 
the normal direction, with unit speed. Since gij(x,t) = f(t)gij(x,0), in this case there is 
only movement in the normal direction, but the velocity X t need not have \X t \ = 1, such a 
requirement puts a restriction on the function /. 

The above shows that, under some assumptions, the idea of thinking of one fixed topological 
manifold and equipping it with a one-parameter of metrics is equivalent to the ideas of chapter |4j 
for when Ai is a hypersurface. 

However, in the topological manifold case with a given time-dependent metric, 

is an interesting example of a parabolic operator on (Ai, (gij(-, t))te[o,T] i n its °wn right, regard- 
less of whether it has any physical meaning. 

For the remainder of this chapter, we will consider one fixed compact, connected, oriented 
and closed topological manifold Ai. Here, closed implies that Ai is without boundary. We 
further assume that Ai is of dimension 1 < n < oo. We will equip Ai with a one-parameter 
family {gij{-, £))te[o,T] of metrics and ask that for each t G [0, T] the map x i-> gij{x, t) is smooth 
and for each x G Ai the map t gij (x, t) is continuous. 

We call (Ai, (■ ,t))te[o,T\ the evolution of Ai and we will be concerned with defining the 



stochastic analogue of (5.2). 

We will see that the new approach to thinking of the evolution of Ai is that the noise will 
be defined on the evolution of Ai and so is much more natural than defining the noise on 
a reference manifold and mapping the noise forward. Further, requiring that t i— > gij(x,t) is 
continuous for every x & Ai, will ultimately allow us to define the notion of a "random metric" 



as presented in section 5.3 



In this chapter we will be using the notation and definitions as presented in chapter [3] 

5.2 A general parabolic stochastic partial differential equation on 
an evolving Riemannian manifold 

As discussed above, we proceed to define the parabolic generalisation of the stochastic analogue 



of (5.2). To this end, fix [/-separable Hilbert space and define 

V t :=H 1 (M,M9t)^), te[0,T} 
H t :=L 2 (M,du(g t y,R), t E [0,T] 
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which can be thought of as the closure of C°°{M) with respect to the norms defined by 



\u\\ Vt 



/ \u\ 2 + |Vm| 2 du(g t ) 
Jm 



Mk : = \l I \ u \ 2 dv(gt) 

M 

respectively. Here, du(g t ) is the Riemannian volume element of M. with respect to the metric 
gij(x,t) and, as in chapter 111 we let \g(x,t)\ denote det (gij(x,t)). 



Let i : U — > H be Hilbert-Schmidt, which by proposition 2.5 always exists and let Ifbea 
[/-valued cylindrical Q-Wiener process with Q = I. 

Since we have that (x,t) i— > gij(x,t) is smooth in x and continuous in t, there exists a 
diffeomorphism 

$ t : (M,g ) — >• (M,g t ) 
for each t E [0, T] and by the smoothness assumptions there exists a', b' > such that 

a' < \D$ t (x)\ 2 < V for every (x,t) E M x [0,T]. (5.4) 



From this, lemma 5.1 below, a change of variable and the chain rule, we see that there exists a 
map F t : Vq — > V t which is bounded, linear and bounded away from 0, uniformly in time and is 
defined in the natural way by 

(F t f)(x)) = fi^t 1 ^)) where x E (M,g t ). 



Observe that F is simply the identity mapping. Indeed, lemma 5.1 and equation (5.4) implies 
that there exists pi,P2, Qi, <?2 > such that 

Pi\\u\\v < \\F t u\\ 2 Vt < qi\\u\\ 2 Vo for every t E [0,T] ^ 
P3||w||fl- ^ \\ F M\ 2 H t < ?2||w||h for every t E [0, T] 

noting that indeed F t makes sense as a map from Hq to H t . 

Lemma 5.1. For t E [0, T] let u : M — > R. Then u E H t if and only if u E H , where in both 
cases M. is equipped with the metric gij(-,t). 

Proof. Since A4 x [0,T] is compact and (x, t) i-> gij(x,t) is continuous it follows that (x,t) i— >■ 
yl^(^ t)| is continuous. Hence, there exists a\, b\ > such that 

ai < v^0M)T < &i for every (x, t) GM x [0, T]. (5.6) 



By Hebey 2000 , for w : Ai — ?■ IR sufficiently smooth 



/ w di/(p t ) = ^2 ( a jV\9(-,t)\w) o ^ 1 dx 

where dx is the Lebesgue volume element on IR n and (Uj,(pj,aj)j & i is a partition of unity 
subordinate to the atlas (Uj,(fj)j e i, that is 

(i) (ctj)j is a smooth partition of unity subordinate to the covering (Uj)j] 
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(ii) (Uj, (fj)j is an atlas of M. and 

(iii) for every j G /, supp(a,) C Uj. 

Note that the atlas of M. is independent of the metric g and so the above charts (fj are 
independent of time. 
Thus, if u G H t then 



f \u\ 2 du(g ) < * / y^tM \u\ 2 du(g ) 



\u\ 2 du(g t ) 



"1 J.M 

< oo. 



Hence 



u\\ 2 Ho < -\\u\\ 2 Ht . (5.7) 



a 



Conversely, if u G i^o then 



|m| 2 dz/((7() < — [ —j= ' = |^| 2 di/(<7t) 

= ~E / ^ 3 ^M\u\ 2 )o^dx 



u\ 2 du(g ) 



= h 

a i Jm 
< oo. 

This completes the proof and shows that 

ll* t <^ll* (5-8) 

□ 

Consider the following parabolic SPDE on (Jvl,g t ) (which can be thought of as the parabolic 
stochastic generalisation of (5.2)) as 

du = Audt + F t idW(t) , 
u(0) = u 

with Gelfand triple V t C H t C V*, where 

Au := div M (cii(Vu)i) — bidiU — cu 
with a,b G C 1 (A / 1; M n ). We assume there exists a, b > such that 

a < ctfc < 6 for every G {!,••• , n} 
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and 



div M (b)<0, beL°°{M). 



We suppose that c G L°°(M) and u G L 2 (Q,F ,f] H ). 

Equation (5.9) is interpreted as an SPDE on (Ai,go) of the following form with Gelfand 
triple VqGHqC V * 

dv = F*AF t vdt + idW(t) 

( 5 - 10 ) 

v(0) = Uq. 



Definition 5.2. Suppose that there exists a solution to (5.10), in the sense of definition 3.16, 
Call the solution v. Then we define the solution to (5.9), u, by 

u{t,u){x) := (F t {v{t,u)))(x) 

where x G (A4,g t ), t G [0, T] and oj G Q. For brevity we shall write u = F t v with the above 
definition in mind. 



The approach above is completely analogous to that of chapter |4| and that (5.10) is com- 
pletely natural, for the reader may verify that F£F t = Ih where Ih is the identity operator 
on Hq. 

The following shows that there is a unique solution to (5.10). 



Theorem 5.3. Let uq G L 2 (Q, Fo, P; Hq). Then there exists a unique solution of (5.10) in the 

sense of definition 3.16 . Moreover, 

E[ sup |Kt)||y <oo. 
te[o,T] 



Consequently, by lemma 5.1 we have 



E[ sup < oo. 

te[o,T] 



Proof. As before, we need to show that F*AF t and i satisfy HI to H4 of assumption 2.2 , and then 
apply theorem 3.17 to see existence and uniqueness. We will see that this is straightforward. 



1. Clearly as A is linear and F t is linear and the composition of linear maps is linear we see 

i is independent of u, v G V we have that 



that F^AF t is linear and so HI is satisfied. 



2. To see H2, let u, v G V . Then as B 
\\B(-,u) - B(-,v)\\ La (u t H ) = and so 

(F t *AF t u - F t *AF t v } u - v) = (AF t (u - v), F t (u - v)) < -a\\V(F t (u - v))\\% 



+ 



I [ div M {b) {F t u - F t v) 2 dv{g t 



+ \\c-\\ L o*\\F t u- F t v\\ 2 Ht 
< g 2 ||c||Loo||« - v\\ 2 Ho 

since diyyn(fe) < 0. So H2 is satisfied with c = 2g 2 ||c||Loo > 0. 
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3. For H3, let v G Vq. Then by the above and using the definition of the VJ-norm 



(F*AF t v,v) < -a(\\F t v\\ 2 Vt - \\F t v\\ 2 Ht ) + \\c\\ L oo \\F t vf Ht . 

Hence 



2(F t *AF t v,v) + \Hi mHo) < cxWvW'^ -c 2 \\v\\ 2 Vo + c 2 

where C4 > exists and ||i||L 2 (t/,.ffo) < C4 as 2 : £/ — > H is Hilbert-Schmidt. So we see that 
H3 is satisfied with a = 2, c\ = 2g 2 (||c||Loo + a) > 0, c 2 = 2p\a > and fit) = c\. 

4. Finally, for H4 let u,v G Vq. Then 

\{F t *AF t u,v)\ < b\\VF t u\\ Ht \\VF t v\\ Ht + \\b\\ L ~ \\ VF t u\\ Ht \\F t v\\ Ht 
+ ||c|| L oo||F t M|| H ||F t i;|| Ht 

< 3max{6, ||6||l<x>, ||c||Loc}||F t u||v t ||F t i;||v t 

< 3gimax{6, \\b\\ L °e, ||c|| ioc }||«||y ||t)||y 

by definition of the V-norm and so ||F t *AFf-u||y o * < C3 1 1 ^ 1 1 v which is H4 where C3 = 
3qi max{6, ||&||l°°) ||c||l°°} > and g(t) = 0. 



Thus, applying theorem 3.17 completes the proof. □ 



Remark 5.4. The uniqueness of a solution to (5.9) is guaranteed up to the choice of map 
F t . Arguably, the F t which we chose is the most natural and there is a natural choice of the 
diffeomorphism $ t which ignores any concept of rotation. 

5.3 A general parabolic stochastic partial differential equation on a 
randomly evolving Riemannian manifold 



Recall section 5.2 where we only asked that (x,t) h-> gij(x, t) is smooth in x and continuous in 
t. This really allows some freedom in the following. 

We still assume that Ai is a compact, connected, oriented and closed topological manifold 
of dimension 1 < n < 00. As an example of a randomly evolving Riemannian manifold, we 
wish to consider the random isotropic evolution of Ai. This means we consider metrics of the 
form 

gij(x,t) := f(t)gij(x,0) xeM 

where gij(x,0) is a given metric such that x 1— > 5^(2,0) is smooth and / is some random 
function, namely the solution of a diffusion equation on K, such that 1 1— > f(t) is almost surely 
continuous. This gives that (x,t) y g^ix^t) is smooth in x and almost-surely continuous in 
t. Equipping Ai with this family of metrics and fixing a realisation of /, gives the random 
isotropic evolution of Ai, whilst retaining the smooth structure of the manifold. 

Recalling that should be positive definite, a suitable choice of / would have that fit) > 
for every t G [0, T], P-a.s. In order for us to mimic the previous section, we want the existence 
of constants a, b > such that 

a < f(u,t) < b for every t G [0,T], P-a.swGSl, 

and so a function / satisfying the above would be preferable. An example of / is given in the 
following. Let (fl, J 7 , P) be a complete probability space and let B(t) be real-valued Brownian 
motion on (Q, J 7 , P). Fix T < 00. Then 
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1. B(0) = P-a.s; 

2. 1 1-)- B(t) is P-a.s continuous; 

3. I? has independent increments. 

Consider the following stochastic differential equation on R, interpreted in the Ito sense 

df(t)=rf(t)dt + af(t)dB(t) 

/(o) = i (5 ' U) 



where r,a e R are constants such that r — ^- < 0. Then by Ito's formula (0ksendal 2003 ) 
one has that 



a 2 



f(t) = exp((r- Y )t + aB(t)). 

Clearly, f(t) > for every t G [0,T] and /(t) < oo for every t G [0, T], P-a.s. We also have 
that t (->■ /(t) is P-a.s continuous and so for each fixed u) G fi we consider the modification of 
B such that 1 f(t) is continuous, so there exists a Q , > such that 

a^j < f(u),t) < b Q for every £ G [0,T]. 

Thus, fix a) G such that t i— >■ /(t) is continuou^J Consider the one-parameter of metrics 
defined by 

sg(x,t) := f{Q,t) gij (x,0) 

where <?y (x, 0) is given and x h> <7y(x, 0) is smooth. We take T < oo sufficiently small so that 
no topology changes occur, such as pinching. Since /(a), ■) is uniformly bounded away from 
and x i->- gij(x, 0) is smooth, we can define the diffeomorphism 

®?:(M,g*)^(M : g?) 

by 

$f (x) := xvTM 

where x G (.M,go)- Note here that .D($^) _1 = l/y/f(u, t) and so an estimate analogous to 
(5.4) holds. Analogous to section 5.2, we define 

V* :=H\M,du(gfy,R), t G [0, T] 
Hf:=L 2 (M,du(gf)-R), t E [0,7] 

and fix U a separable Hilbert space, letting i : Z7 —> Hq be Hilbert-Schmidt. Let W be a 
[/-valued cylindrical Q- Wiener process with Q = I. We define the natural map between V$ 
and by 

Jf:Vf— >Vf (Ff M )(x):= M (x/v^t)) 

where x G (.M, g^). This defines a natural map between if^ and Hf and an analogous inequality 
to (5.5) holds, with the Pi,qi, (i = 1,2), dependent on Co. 

We now equip M. with this one-parameter family of metrics (gfj(-, £))te[o,T] an d consider the 
following general parabolic SPDE on (Ai,gf) 

du = A°udt + FfidW(t) 
u{0) = w 



4 By which we mean we consider the continuous version of the Brownian motion (which exists by Kolmogrov's 
continuity theorem (0ksendal 2003 , theorem 2.2.3)) and fix a realisation. 
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with Gelfand triple Vf C Hf C Vf *, where 



A°u := div^ ( (a i (V £i 'u) l -) - hdfu - cu 



with a,b G C 1 (A^; IR n ) and there exists a, b > such that 

a < ak <b for every G {1, 

and 



di<,(6)<0, 6gL°°(7W). 

We suppose that c G L°°(.M) and uq G L 2 (f2, J^, P; -f^o )• Here we make the implicit assumption 
that for any realisation uj, div^(6) < 0. 

As before, equation 5.12 is interpreted as an SPDE on (M,g^) of the following form with 
Gelfand triple V? C H$ C Vg* 



dv 
v(0) 



F?*AF?u dt + idW(t) 

Uq 



(5.13) 



and the solution to (5.12) u, is defined as (cf definition 5.2) 

u(t,u)(x) := (F?(v(t,u)))(x), 
where t G [0,T],x G (M,gf) and w G f2. 

Theorem 5.5. Let uq G L 2 (fi, .Fq, P; -f^o)- 27ien there exists a unique solution (in the sense of 



definition 3.16) to (5.13) and 



E[ sup IKOIiy <oo. 
te[o,T] 



Consequently, by lemma 5.1 we have 



E[ sup || M (t)|| 2 ,] <oo. 
te[o,T] 



Proof. The proof is completely analogous to that of theorem |5.3| noting that in our case cu is 
fixed and all the properties we exploited in the proof of theorem 5J3 also hold here, although 
some of the bounds are dependent on the underlying probability space Q. To save needless 
repetition the reader is directed to the proof of theorem 5.3| □ 

Remark 5.6. We chose random isotropic evolution of Ai to illustrate that we need only that 
t i — y gij(-,u,t) is continuous for each fixed Q G Q where we consider the continuous version of 
the Brownian motion such that t i— >■ fit) is continuous, fixing a realisation. Of course, we still 
need x h-> gij(x,u, •) to be smooth. Since the initial metric is chosen sufficiently nice, this also 
holds. 

However, the above results can be easily seen to hold for general random metrics gij(x,u,t) 
where u is fixed, is smooth in x and continuous in t. Above all, it is important to ensure the 
existence of a diffeomorphism <E>^ : {Ai^g^) —> (Ai,gf) and arguably random isotropic evolution 
yields the simplest example where one can write down 

With the assumptions outlined above, we see that no topology changes as long as T is chosen 
sufficiently. However, it is relatively easy to come up with examples where topology changes can 
occur. 
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Example 5.7 (An example of a topology change). Let n = 2 and consider the level set of the 
function 

Figure^ graphs the level sets of f(x,y) = c for c = 0.3,0.4,0.5 and 0.6. We see that as c 
decreases to 0.5, the smooth curve becomes pinched. As c decreases further, two distinct curves 
are produced and are clearly not diffeomorphic to the level curve at c = 0.6. This shows that a 
topology change has occured. 

Although the above example is deterministic, one can imagine a random isotropic pertur- 
bation of the c = 0.6 level set where the perturbation is bounded above and below, but with 
sufficiently high values as to cause pinching like that of the c = 0.5 level set. Such a scenario 
was not considered in this chapter. 




Figure 1: Graph of the level sets of / for c = 0.3, 0.4, 0.5 and 0.6. 
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6 Further research 



Although we have formulated what it means to have a SPDE on a moving hypersurface and 
then looked at SPDEs on an evolving manifold, there are still many avenues of research to 
pursue for the future. 

The approach that has been used is the so-called variational approach. This approach 



is not widely used, mainly due to the constraints of HI to H4 of assumption 2.2 which are 
perhaps too restrictive in certain circumstances. For example, when one takes A := A — f 
where f(s) = s 3 — s we see that A no longer satisfies H2 of assumption 2.2 for / is not globally 
Lipschitz, nor monotone. Indeed, taking A := AA leads to the Cahn-Hilliard-Cook equation 

amongst others). 

where 



(Da Prato and Debussche 1996 , Kovacs et al. 2011 



1992 



Perhaps a more natural approach would be that of Da Prato and Zabczyk 
weak solutions (in the sense of weak solutions to PDE) are considered. This is the approach 
used by Gyongy [1993 for formulating SPDEs on differentiable manifolds. Hence, our work is a 
generalisation of this for the existence and uniqueness theory for the variational approach. An 



ideal research avenue would be to repeat the above theory, but for the Da Prato and Zabczyk 



1992 approach. 



The noise that was considered was a ^/-valued cylindrical Wiener-process. There are many 
more examples of "noise" to be considered. For example, white noise, which is constructed on 
the Schwartz space of tempered distributions as in Holden et al. 2009 could be considered. 



A problem here is that it is not obvious what the generalisation of the Schwartz space of 
functions on M n is for an arbitrary Riemannian manifold. Indeed, such a generalisation exists if 
the manifold is a Nash manifold, Aizenbud and Gourevitch 2007, 2010 . Another type of noise 



that is becoming popular in the stochastic analysis community is a class of Levy processes. 



Peszat and Zabczyk 2007 



Here the Da Prato and Zabczyk 1992 method is extremely useful, as this is the method used 
in 



Here, the infinite dimensional Levy process has already been 
constructed and so what remains is to formulate, analogously to this paper, SPDEs on evolving 
manifolds which are driven by a Levy process and to formulate the variational approach to the 
analysis of SPDEs. 

The above are a few possibilities for generalisation for when the metric of the manifold is 
evolving deterministically, or indeed randomly. However, the ultimate goal is to consider when 
the metric evolves depending on the solution of the SPDE on the manifold. This notion can be 



found Neilson et al. 2010 for a coupled system of SPDEs where the hypersurface evolution is 



coupled to the solution of the SPDE. 

This is an extremely challenging mathematical problem. To attack it, one must first be able 
to understand what the actual problem is mathematically. For example, if we fix to G Q and 
consider the metric given by gij(u", •, t) then we are in the position of section 5.3, however, now 



-==L===JL (g i i{u",x,t)y/\g{u«',x,t)\^- 
^\g(u^,x,t)\dxi V Ox, 



in local coordinates, which is a nonlinear differential operator on Ai. We are now considering 
the equation (dropping the superscript u) 



du = AJ^udf + zdW(t) 
u(0) = u 



(6.1) 



where i is some suitable Hilbert-Schmidt operator. Since u is fixed, the above is in fact a 
nonlinear PDE on H, whatever H should be. Naively, we would set H := L 2 (Ai, dv(gt(u)); M) 
so now the Lebesgue and Sobolev spaces also depend on the solution u of ( 6.1[ ). 
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Suppose we have isotropic evolution, so the metric is given by gij(u, x, t) = f(u(u, t))gij(x, 0) 
where / : H — > R is such that there exists a, b > such that a < f(u) < b for every 
u G H . Then, if the initial metric is sufficiently nice, we have that a' < ^J\g{u)\ < b' for 
suitable a', b' > 0. So, assuming no topology changes occur, we have to establish existence and 
uniqueness theory for the non-linear operator in this space, noting that we do not have 
that H C L 2 (A4, \/\g(u)\ di/(<7o)> R) 5 which would be helpful in the attaining of estimates for 
HI to H4. Note also that the definition of H is circular and therefore we should be looking 
for another space of solution, or even a different notion of solution such as viscosity solutions 
(Evans 1998| ). Since (6.1) is purely deterministic, a first approach to this problem would be 



looking at the deterministic analogue. Unfortunately this theory does not exist. 

The reader will note that throughout this paper, all the functions are real valued. In 



Funaki 1992 SPDEs with values in a Riemannian manifold are mentioned. This is really 



a generalisation of the pioneering work of El worthy 1982 , where SDEs on manifolds were 



considered. A natural extension to this paper would be developing the above theory for SPDEs 
whose solution is a function with values in the manifold. The approach of chapter [5] will be useful 
here but in order to formulate the SPDE analogously to the above, we would have to consider 
C°°(S x TAA) valued Wiener process, where S is the unit circle and TM. is the tangent bundle 
( |Funaki| p92] ). However, C°°(S x TM) is not a Hilbert space and the equations presented 



are in the Stratonovich form (not the ltd form as in this paper), so the theory of chapter M 
would have to be adapted to this space, with a suitable topology. This would then yield a 
theory of the variational approach to SPDEs whose solutions are functions taking values in 
the manifold, where the metric on the manifold is given by a one-parameter family of metrics. 
Further generalisation could be found by asking that this family is random as in chapter [5j 

In this paper, only specific examples of operators were considered, until chapter |5j This 
was really to build up intuition as to what spaces one should set the equations in. Using 
the theory of chapter |5j an extension to this project would be considering general parabolic 
operators, with certain non-linearities, on M where M has a one-parameter family of random- 
metrics associated to it. It is fairly obvious how to give a time-dependent generalisation of the 
assumptions HI to HA to give a general variational theory of SPDEs on evolving Riemannian 
manifolds. 

In all the above extensions, only existence and uniqueness properties have been discussed. 
What would be interesting to look at is long time behaviour of solutions. For technical reasons, 
we have only considered a fixed finite time T, and assumed that the topology does not change 
over [0, T]. To extend, topology changes could be considered and long time behaviour of the 
solution on the manifold could also be considered. In light of the above discussion, if one is able 
to say something about the solution of a SPDE whose solution is a function taking values in 
the manifold where the solution is coupled to the evolution of the metric, one could perhaps say 
something interesting about the long term behaviour of the manifold and the solution. This 
would be truly remarkable, because then one would have (hopefully) some convergence results 
of manifolds and the solution. 

In conclusion, there is much scope for developing the mathematical ideas presented in this 
paper for future research. 
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